Dynamical density correlation function of ID Mott insulators in a magnetic field 
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We consider the one dimensional (ID) extended Hubbard model at half filling in the presence of a 
magnetic field. Using field theory techniques we calculate the dynamical density-density correlation 
function Xnn(u),q) in the low-energy limit. When excitons are formed, a singularity appears in 
X„„(u,(/) at a particular energy and momentum transfer. 



I. INTRODUCTION 

Quasi ID Mott insulatorsEl display unusual phenomena 
like spin-charge separation and dynamical generation of 
a spectral gap and have therefore attracted much atten- 
tion in recent years. At present the best realizations of 
ID Mott insulators are found in anisotropic antiferromag- 
nets like SrCuC>2 or Sr2Cu03. The dynamics of the lat- 
ter compound lias been studied both by angular-resolved 
photoemissioncl and by electron energy-loss spectroscopy 
(EELS)H. EELS measures the dynamical density-density 
correlation function. Theoretical descriptions of ID itin- 
erant antiferromagnets are based on the half-filled Hub- 
bard model with a large on-site repulsion U (compared 
to the hopping matrix element t), which ensures that the 
single-particle Mott gap is large. In the case of S^CuOa 
it is believed that U ~ 8t and that density-density in- 
teractions batween neighboring sites need to be taken 
into accountBu. Given that U is "much" larger than t, 
strong-coupling expansions around U = 00 are an ap- 
propriate starting point. They have been used to deter- 
mine dynamical correlation functionscra and to success- 
fully model the EELS data for SrCuOjl 

There exist other materials believed to be quasi-lD 
Mott insulators (e.g. the Bechgaard saltsH), in which the 
Mott gap is small compared to t. The "weak-coupling" 
regime U < 2t in which the Mott gap becomes small is 
manifestly beyond the range of applicability of strong- 
coupling expansions. The very existence of a-i-sap pre- 
cludes the application of conformal field theoryQcL How- 
ever, this regime is accessible by an approach based on 
exact field theory methodHtM In Refs [T3iiMlp| this 



method has been used to determine the optical conduc- 
tivity and the single-particle Green's function for ID 
Mott insulators in the weak-coupling regime. Here we 
employ this approach to calculate the dynamical density- 
density correlation function. Motivated by suggestions 
that a magnetic field may generate X-ray edge like thresh- 
old singularities!!^ in the density-density response, we 
take into consideration the effects of a magnetic field. 
The outline of this paper is as follows: in section II we 
present the field theory description for the low energy 
degrees of freedom. In section [Hi] we determine the dy- 
namical density-density correlation function Xnn(w '» (?) • 
In section [V we consider an extended Hubbard model 



density-density interactions. Here excitons are formed 
and we determine their contributions to Xnn{u,q). Wc 
summarize our results in section [y| 



II. FIELD THEORY DESCRIPTION 

The extended Hubbard model in a magnetic field is 
described by the Hamiltonian 



3=1 t 



where c; i<T are fermionic annihilation operators of spin f 
(cr = 1) and J. (cr = -1), n Lcr = cj )0 .c/ >er and ni = n; iT + 
ni,i- For h = the ground state has spin projection S z = 
whereas for very large fields h > h c it is fully polarized. 
We constrain our analysis to the case where the ground 
state is partially magnetized, which corresponds to fields 
< h < h c . A description of the low-energy degrees of 
freedom of ([!]) for weak interactions < 2V2 r< 2Vi < 
U <C t is then obtained by standard techniquesO. In the 
presence of the field h there are four Fermi points ±kp t(T 
with Jcf,i + = 71 '/ /fl 0j where ao is the lattice spacing. 
Taking into account only modes in the vicinity of k,F,a 
we expand 

<*,* — » >/5o b lkF -° x + e- ik "'" x L a (x)] , (2) 

where x = IdQ. The resulting fermionic field theory can 
be bosonized with the result 



£ s = 



1 

167T 



v s (d x $ s ) 2 - -(ft$ s ) 5 

V. 



£ c = ^- \v c (d x $ c ) 2 - -(d t $ c ) 2 

f07T V c 



(3) 

Acos(/3 c $ c ).(4) 



The spin sector is a free bosonic theory whereas the 
charge sector is described by the integrable Sine-Gordon 
model (SGM)I3. Fermionic operators are expressed in 
terms of the canonical charge and spin bose fields & c ,s 
and their respective dual fields 



with sufficiently strong nearest and next-nearest neighbor 



-1 f x 

Q c>s (t,x) = / dy d t $c,s(t, 

Vc,s J — 00 



y) 



(5) 
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by 

L a = i la e i(^-* e =) el^.*.-*®-), 

= ^ e -i(^*«+7fe e «) e -i'0»-«-+7fce.). (6) 

Here r? CT = 77J are Klein factors that fulfill {fyo-, 7 ?-;-} = 
25 CT)T . The spin and charge velocities v c , s and the param- 
eters /3 CjS depend on t, h, U and V. 

A. Hubbard model (Vi,2 = 0) 

For the Hubbard model /3 C:S and v CjS can be calculated 
exactly from the Bethe Ansatz solution. The "77-nairing" 
SU(2) symmetry of the half-filled Hubbard modeE3 fixes 
(3 C = 1, whereas (3 s (h, U) iaxibtained from the solution of 
a linear integral equatioroia 



V2Z(A) 



Z{\) = 1+ / dfia 2 (X-n) Z(n) 



(7) 



where 27022(2:) = U[x 2 + U 2 /A\~ x . The integration 
boundary A is determined by the condition e(A) = 0, 
where 



e(A) = h - 4RcV 1 — (A — iU/4) 2 + U 

+ [ du. 02(A-/*) e(/*) • (8) 
J —A 

In the limit U/h — > f3 s can be calculated analytically 
1.4 
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FIG. 1: /3s for the Hubbard model as a function of /i for 
different values of U. 



by means of a Wiener-Hopf analysis 

U 



0s =S 1 + 



167rf cos(7rA / J) 



(9) 



where M. is the magnetization that is calculated from the 
solution to an integral equation similar to (H) (see Ref. 
|h. We note that j3 s varies between 1 for h = and \/2 
for /i — > ft, c . The behavior of /3 S as a function of /i for 
different values of U is shown in Fig|l[ We see that for 
small values of U (3 S remains very close to 1 up to large 
fields very close to h c . 

The spin and charge velocities can be calculated in a 
similar way. Field theory is exact in the scaling limit, 
which has been constructed in the absence of a mag- 
netic field in Ref-gO using exact results for the half-filled 
Hubbard modeEST The h > case can be mapped 
onto the attractive Hubbard model below half-filling by 
means of the particle-hole transformation for spin down 



■1)3 



The scaling limit for the latter model 



has been found by Woynarovich and ForgacsE3 and is ob- 
tained by taking t — > 00, U — ► while keeping 

y/Ut cos 3 (irM)e- wtcos ^ M ^ 2U = fixed. (10) 

In this limit (3 S = 1 and the low-energy effective field the- 
ory is SU (2) x SU(2) symmetric. However, on the level of 
the large-distance asymptotics of correlation functions of 
the underlying Hubbard model this enhanced symmetry 
is broken down to SU (2) x U(l) by the oscillating factors 
in (Q). For example, the leading asymptotical behavior 
of the spin-spin correlation functions is 



(S x (x)S x (0)) = (S v (x)S y (0)} 



B(-l) x / a « 



(S*(x)S*(0)) 



Acos2kp a x 



(11) 



and the spjp=jS'C/(2) symmetry is broken. Previous 
experienceEj'ta suggests that field theory gives a good 
description of the lattice model in an extended vicinity 
of the scaling limit, provided the gap is small compared 
to t. This is the case as long as U < 2t. We will apply 
field theory in the regime defined by this criterion and 
therefore allow (3 S to be different from 1. 



B. Extended Hubbard model 

For Vi,2 ^ no exact results for the lattice model (|l|) 
are available. However, the low-energy degrees of free- 
dom of (|l|) in the regime < V2 < V\ < ^ <t are 
still described by (|)-(|), but now with C < ill Re- 
cently it was shown in Ref.|l4| that for sufficiently large 
values of U, V\ , V2 it is possible to reach the attractive 
regime of the SGM (3 C < l/y/2, in which excitonic holon- 
antiholon bound states form. We will first consider the 
range l/v2 < j3 c < 1, where no excitons eptist and field 
theory results for the optical conductivityliia have been 
found to be in good agreement with dynamical density 
matrix renorrftalization gro up computations for the lat- 
tice model (0)cil. In section IV we extend the analysis to 
the regime (3 C < 1/ \f2. 
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C. Density Operator 



The density operator is expressed in terms of the spin 
and charge bosonic fields asLj 

n(x,t) = n (a, t) + y~] n 2 k F ,, (x, t) , 

n {x,t) = A d x $ c , 
n 2kp Jx,t) = A' e 2ik ^ x sin(^$ c ) e^*= . (12) 

Here A and A' are numerical constants. For a less-than 
half-filled band there is an additional contributiontZl 



nu(x, t) = Ay cos (2[/cf,t + + /?c$ c ) 



(13) 



which is obtained by integrating out the high-energy 
degrees of freedom in the path integral representation 
for the density-density correlation function of the lattice 
model. The operator ( |l3"| ) corresponds to scattering pro- 
cesses involving two particles and two holes with momen- 
tum transfer 2(kp^ + fcjrjj = 2n/aQ. As a result Ajj is 
proportional to U/t and thus is small. At half-filling we 
have Ajj = 0. This can be established by considering the 
following discrete symmetry of the lattice model 



(14) 



The lattice density operator transforms as rij — > 2 — n 3 . 
In the field theory this symmetry corresponds to inverting 
the signs of the charge boson and its dual field <& c — > 
— $ c , 6 C — * — 6 C - The (normal ordered) density operator 
n(x,t) must transform to — n(x,t) under this change of 
sign and this implies that Ajj — 0. 



III. DENSITY CORRELATIONS FOR §</?;?< 1 

The density-density correlation function is given by 

G nn (x, t) = G° nn + 2_] G 2 nn' a i (15) 



where 



G 



nn 



\n (x,t) n (0,0)|0) , 
\n2k Fa {x,t) n 2 fc F „„(0,0)|0) . 



(16) 



We note that there are no "mixed terms" as can be shown 
by exploiting the transformation properties under charge 
conjugation. Due to spin-charge separation the corre- 
lation functions in (|l^) factorize into spin and charge 
pieces. In the spin sector (^) we are dealing with a sim- 
ple Gaussian model and elementary considerations give 



,<0|e 



>|0).= [x 2 - (V s t + It) 2 ] 



(17) 



where |0) s denotes the vacuum in the spin sector and 

d = 0112. (18) 



In the Hubbard model the exponent d varies between i 
for zero field and 1 for h — ► h c . In order to determine 
the charge part of the correlators ( |l6| ) we make use of 
the integrability of the SGM (^) describing the charge 
sector. In the range of (3 C considered here the spectrum 
of the SGM consists of scattering states of solitons and 
antisolitons, which are particles of mass M, charge Q = 
±e and relativistic dispersion e(p) = \Jp 2 + M 2 . In the 
Hubbard model they correspond to holons and antiholons 
respectively. It is convenient to parametrize energy and 
momentum in terms of the rapidity variable 9 



M 

p = — sinh 9, e — M cosh 6. 

v c 



(19) 



We introduce an index e — ± for solitons and antisoli- 
tons. Then a scattering state of n solitons/antisolitons 
with rapidities and internal indices {£&} is denoted 

by \9 n . . . 9i) En £l . In the spectral representation of this 

basis of (anti)soliton scattering states we may express the 
two-point function of an operator O in the charge sector 



: (0|O(x,i)O t (0)|0) c = EE 



n=0 e 



d9 1 . . . d9 n 



x exp 



3 = 1 



(2n) n n\ 
<O|O(O)|^...0 1 ) e „... £l | 2 .(2O) 



Here pj and ej are given by (|l9|), and the form factors 
c (0\O(0)\9 n . . ■9i) en ... £l \ can b£,-calculated by exploiting 
the integrability of the SGMEJ'EJ. As a consequence of 
the transformation properties under charge conjugation 
of the operators appearing in (|T^), only intermediate 
states with an even number of particles will contribute to 
(pp|). In order to obtain an accurate result for the large- 
distance asymptotics it is sufficient to take into account 
intermediate states with only, a small number of parti- 
cles in„ihe spectral sum (20)123. For the case at hand we 



havi 



|sin(^)|^ 2 ) ±T | 2 



|/(2^)sinh0 + | 2 ; (21) 
Z 1 \f(29^)\ 2 , (22) 



where 9± 

m 

F(9) 



9\ ± 6*2 )/2, Z\ is a known constantS and 
F{6) 



cosh 



( 9+ir L \ ' 
V 2 « J 



sinh- exp 
2 



dk sin 2 (|[f +i])sinh(i^fc) 



sinh 



(f) 



cosh (■§ ) sinh k 



(23) 



EELS measures the imaginary part of the retarded dy- 
namical density-density correlation function 



Im< i 



dx I dt e 
la 



itjt—ikx 



Grin \X, t) 
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(24) 



The threshold is obtained by minimizing the energy at 
fixed total momentum with respect to p, q\ , q2 



We evaluate Xnn( w >^) i n the vicinity of the low-energy 
modes at k = 0, 2kp i<7 by Fourier transforming the large- 
distance asymptotics of the density-density correlation 
function. The latter is obtained by carrying out the form 
factor expansion ( p0| ) in the charge sector and multiply- 
ing it by the spin-piece ( |l7j ) in the case of the 2kp a re- 
sponse. 



A. Small k behavior 

In the vicinity of k = the dynamical density response 
is dominated by the contribution from G^ n (x,t). This 
contribution does not involve the spin sector and straight- 
forward calculations give 



Xnn(w,g) = A 2 8M' 



M 2 |/o(0o)f 



AM* 



9(s 2 -4M 2 ), 



(25) 



with s 2 {w,q) — uj 1 — (v c q) 2 and 9 — 2arcosh(s/2M). 
Above the threshold at id = \J v 2 q 2 + AM 2 , Xnn(u,q) 
increases from zero in a universal square root fashion. 
This is due to the momentum dependence of the form 
factors. 

From the Heisenberg equations of motions for the lat- 
tice density operator one can derive a relation between 
Xnn and the optical conductivity a(tu) 



Re a(u>) = lim — q) 

q^O q l 



(26) 



The optical conductivity has been calculated inE3 and 
agrees with (J2q ) and (£5|). We note that a contribution 
of the type (|13|) to the density operator would violate the 
relation (^6|)_-which is an independent argument showing 
that A v = 0Ei 



B. Behavior around k — 2/cf,ct 

In the vicinity of k = 2kp. a the dynamical density re- 
sponse is dominated by the contribution from G^n'" (x, t) 
and involves both the spin and the charge sector. The 
threshold can be determined by considering the lowest 
intermediate state that couples to the density operator 
at k = 2kp ai which is a scattering state of one soliton, 
one antisoliton and one spinon. The total momentum 
and energy of this state are 



P = p + qi + q2 , 



E = v s \p\+^2^l M 2 + v 2 q 

i=i 



2„2 



(27) 



mm 

91.92 



V s \P~qi-q2\+J2y M2 + V c1* 
if \P\ < Q 



\v s \P\ + 2MVl- a 2 if|P|>Q' 



where 



v s ^ 2Mv s 
v c v c\/v 2 -v 2 



(29) 



The behavior is rqiiite similar to what is found for the 
spectral functionlia. 



1. Equal velocities v a — v c — v 

For v s — v c = v one can obtain the following represen- 
tation 



Xnn(^,2fep> + q) 



Y 2 {l-d)Z 1 A" 



2<2-l 



x ImF 1 - d, 1 - d, 1 



tt(2v) 

2 2 2 

uj z — v q 
' c 2 (6) 



de 



\hm\ 2 

c (0)2-2d 

(30) 



where c{9) = 2Mcosh(0). The imaginary part of 
the hypergeometric function vanishes unless c(9) < 
\/ oj 2 — v 2 q 2 . This implies that the response function is 
nonzero only if uj 2 > v 2 q 2 + 4Al 2 , in agreement with (|2q) . 
Just above the threshold (0 < s/2M <§; 1) we may use the 
transformation formulas for hypergeometric functions to 
obtain (s 2 = uj 2 — v 2 q 2 ) 



ImF 1 - d, 1 - d, 1, 



c 2 (9) 



xF d,d,2d, 1 



c 2 (9) J \c 2 (9) 



x e(s 2 -c 2 (9)), 



r(l-2d) sin(7r2d) 
T 2 (l-d) 

2d-l 



(31) 



where 0(x) is the Heaviside function. The remaining 
^-integral in (|3^) is therefore over a very small interval 
[— arccosh(s/2M), arccosh(s/2M)] and can be taken by 
Taylor-expanding the integrand. The leading contribu- 
tion to the behavior just above the threshold is 



Xnn(^,2fc Fj(T + q) cx 



s-2M\ 2 



M 
s - 2M 



V'2d 



M 



0. 



(32) 



Here we have used that fl 2 > 1/2, in which case we have 
/i(20) oc0 , for 0-> 0. (33) 
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At the Luther-Emery point 0% = ic3 there is a different 
power law increase in ( |3^ ) (the exponent is 2d — |). As 
— * ~ from above the region is which (|3^) holds shrinks 
to zero. The important result here is that Xnn vanishes 
as the threshold is approached from above. There are no 
threshold singularities! The behavior for large frequen- 
cies u> 3> \J v 2 q 2 + AM 2 (but necessarily w < t for field 
theory to apply) is 

X„„(w, 2fc F)ff + q) oc fl M - 2 +^ , s » M. (34) 
2. Different velocities v s 7^ v c 



spectral function, that are associated with v s and v c 
respectivelyatJ. A natural question is whether an anal- 
ogous feature exists in the density-density response. We 
have analyzed (j3f]) for several sets of parameters f3 c , (3 S , 
a and found that in all cases x„ n (w, 2kp a + q) is rather 
featureless. There are no singularities or peaks that can 
be associated with v c and v s separately. We also do 
not find any threshold singularities as the magnetic field 
is increased. In Fig. ^ we plot x„„(w, 2kp ia + q) for 
j3 c = (3 S — 1 and a — 0.851, which corresponds to the 
half-filled Hubbard model in zero magnetic field at U = 1 
(v c = 2.15tao, u s = 1.83tao). At uj ~ 8M one can just 
see that the threshold is below the curve -J 'AM 2 + v 2 q 2 . 



In the case of different spin and charge velocities ^f- 
a < 1 one may represent %„„ as 



Xnn(u,2k F>(T + q) 



Z X A! 1 



r 2 (d)(2v s ) 2d - 1 



(35) 



where 



tt = uj - v s q - 2Mcosh(0_) [cosh(0 + ) - a sinh(0+)] , 
fi' = uj + v s q-2Mcosh(6-)[cosh(d + )+asxnh(9 + )]. 

(36) 



One easily checks that (|35|) leads to a threshold described 
by (28). The remaining integrals in ( pS5| ) are evaluated 
numerically. 



cr 4 
+ 



3 



2 - 




4 6 

co/M 
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IV. EXCITONS: ^ < § 



In the regime f3 2 < ^ soliton and antisoliton can form 
excitonic bound states which for the SGM are known as 
breathers. There are 



N = 



l-(3 2 



(37) 



different types of excitons, where [x] denotes the integer 

part of x. We denote the different excitons by e\, e^, 

The exciton gaps are given by 



M n = 2M sin(mr£/2) , n=l, 



(38) 



It follows from the transformation properties under <i> c — > 
— $ c that only the "odd" excitons e2n+i couple to the 
density operator. For simplicity we take f3 c > 1/2 from 
now, in which case we only have to consider the first ex- 
citon e\. In this regime of (3 C the leading contributions of 
the spectral sum to the dynamical density-density corre- 
lator are given by 



Xnn(w, k) = X° X n(w, k) + Xnn( W > fc ) 



(39) 



where x™(^,^) and x«n(w,fc) denote the contributions 
from intermediate states with one exciton and many 
spinons and one soliton, one antisoliton and may spinons 
respectively. We have already calculated Xnn( w 'fy 
above. For small k it is given by (pjj) and for k w kp, a 
by (|35|). The exciton contribution can be calculated by 
the same methodEi and we now present the results. 



FIG. 2: X nn(Lo,2k F ,« + q) for fl. = fa = 1, a = 0.851 
and several different values of v c q. The curves have been 
offset. The dashed line is the threshold when spinons carry 
zero momentum + via . 

A remarkable feature of the half-filled Mott insula- 
tor is the presence of two dispersing features in the 



A. Behavior around k — 

Here the exciton is visible as a sharp (^-function peak 
at an energy below the soliton-antisoliton scattering con- 
tinuum 

X$f„>, q) = A 2 g ^6(u - Jv 2 q 2 + M 2 ) , (40) 

UJ v 
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where 



.90 



2 £2 




2 sin | ^ | c.\j) 




(41) 

The result ([l0|) is again related by the equations of mo- 
tion ( p6[ ) to the corresponding contribution to the op- 
tical conductivitylij. The dynamical density suscepti- 
bility @, for w c g = 0.2M, is plotted in Fig. [| For 
f3 c < l/\/2 the breather contribution ( |4C| ) appears and 
the spectral weight is gradually transferred from the 
soliton-antisoliton continuum to the coherent peak. 
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FIG. 3: Xnn(w,q = 0.2M/v c ) for different values of /3 C . 
We have broadened the delta function by convoluting it 
with a Lorentzian in order to exhibit the transfer of spectral 
weight from the soliton-antisoliton continuum to the coherent 
breather peak. 



The exciton contribution to the dynamical density 
density correlation function is given by 



xE£(w, 2k F , a + q) 

K 1 



T 2 (d)(2v s ) 2d - 1 



2cos«/2) 



/oo 
^(ss') d_1 e(E) e(s'), 
-oo 



(44) 



where 



E = oj — v s q — Mi [cosh(0) — a sinh(0)] , 
E' = uj + v s q - Mi [cosh(0) + a sinh(0)] . (45) 

One readily deduces by inspection of equations ( f45|) and 
( Eij ) that for hxed q the exciton contribution to x exhibits 
a cusp at a frequency 



Wcusp = \l Ml + v 2 q 2 



(46) 



If we fix the momentum transfer to be q = Q' , this cusp 
turns into a singularity. In order to exhibit these interest- 
ing features we plot x„„(<jj, 2kp,a + q) for £ = 0.6, a = 0.8 
and different values of q in FigjiJ 
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B. Behavior around q — 2/cf,<t 

In the vicinity of q = 2kp^ the exciton contributes to 
the dynamical density-density correlation function via an 
exciton-spinon scattering continuum with threshold 



Ethres — 



mm 



v s \P~q\ 



M 2 + v 2 q 2 



| y/M't + v 2 P 2 if \P\ < Q' 



\v s \P\ + Mi VI- a 2 if \P\>Q 



, , (42) 



where 



Q' = 



Miv 



1V S 



(43) 



FIG. 4: x»«(^> 2fc F , CT +g) for a = 0.8, f3 s = 1 and £ = 0.6. For 
ui ~ Wcusp one can clearly see the cusp due to the breather 



contribution. As v c q approaches v c Q' 
turns into a singularity. 



2.157M the cusp 



V. SUMMARY 

We have studied the dynamical density-density re- 
sponse of half-filled ID Mott insulators for the case where 
the Mott gap is small compared to the hopping matrix el- 
ement t. We have allowed for the presence of a magnetic 
field that partially magnetizes the ground state. Unlike 
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the spectral function, the density-density response func- 
tion does not exhibit prominent, dispersing features as- 
sociated with the spin and charge degrees of freedom re- 
spectively. Due to the momentum dependence of the ma- 
trix elements in the gapped charge sector, Xnn(w, k) tends 
to zero as the threshold is approached from above: irre- 
spective of the magnitude of the applied magnetic field 
there are no threshold singularities. 

In the parameter region of the extended Hubbard 
model where excitons are formed, the dynamical density- 
density response exhibits a cusp at some specific value of 
energy transfer for momentum transfers k close to 2kp,a- 



This cusp turns into a singularity for one particular value 
of k. These features should be experimentally observable 
in small-gap quasi-lD Mott insulators. 
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